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I The absolute visual magnitude as function of the observed color (B-V) , also 

> named Hertzsprung-Russell diagram can be described through five equations; 

^ that in presence of calibrated stars means eight constants. The developed 

l/^ framework allows to deduce the remaining physical parameters that are mass 

, radius and luminosity. This new technique is applied to the first 10 PC , 
^ the first 50 pc , the Hyades and to the determination of the distance of 

a cluster. The case of the white dwarfs is analyzed assuming the absence 
O of calibrated data: our equation produces a smaller in respect to the 

K*- standard color-magnitude calibration when applied to the Villanova Catalog 

^ of Spectroscopically Identified White Dwarfs. The theoretical basis of the 

H formulas for the colors and the bolomctric correction of the stars arc clarified 

through a Taylor expansion in the temperature of the Planck distribution, 
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1 Introduction 



The diagrams in absolute visual magnitude versus spectral type for the stars 
started with Hertzsprung (1905); Rosenberg (1911); Hertzsprung (1911) . 
The original Russell version can be found in Russell (1914b, a,c), in the follow- 
ing H-R diagram. The common explanation is through the stellar evolution, 
see for example chapter VII in Chandrasekhar (1967). Actually the presence 
of uncertainties in the stellar evolution makes the comparison between theory 
and observations an open field of research , see Maeder & Renzini (1984); 
Madore (1985); Renzini & Fusi Pecci (1988); Chiosi et al. (1992); Bedding 
et al. (1998). Modern application of the H-R diagram can be found in de 
Bruijne et al. (2001) applied to the Hyades when the parallaxes are pro- 
vided by Hipparcos, and in Al-Wardat (2007) apphed to the binary systems 
COU1289 and COU1291. 

The Vogt theorem , see Vogt (1926), states that 

Theorem 1 The structure of a star is determined by it 's mass and it 's 
chemical composition. 

Another approach is the parametrization of physical quantities such as 
absolute magnitude, mass , luminosity and radius as a function of the tem- 
perature , see for example Cox (2000) for Morgan and Keenan classification 
, in the following MK, see Morgan & Keenan (1973). The temperature is not 
an observable quantity and therefore the parametrization of the observable 
and not observable quantities of the stars as a function of the observable 
color is an open problem in astronomy. 

Conjecture 1 The absolute visual magnitude My is a function, F, of 
the selected color 

Mv^F{ci, ,cs,{B-V)) . 

The eight constants are different for each MK class. 

In order to give an analytical expression to Conjecture 1 we first analyze 
the case in which we are in presence of calibrated physical parameters for 
stars of the various MK spectral types , see Section 2 and then the case of 
absence of calibration tables, see Section 3 . Different astrophysical envi- 
ronments such as the first 10 pc and 50 pc , the open clusters and distance 
determination of the open clusters are presented in Section 4. The theoreti- 
cal dependence by the temperature for colors and bolometric corrections are 
analyzed in Section 5. 
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2 Presence of calibrated physical parameters 



The My ,visual magnitude, against {B — V) can be found starting from five 
equations , four of them were aheady described in Zaninctti (2005). When the 
numerical value of the symbols is omitted the interested reader is demanded 
to Zaninetti (2005). The luminosity of the star is 

logio(7^) = O.4(Mb„i,0-Mboi) , (1) 

where M^^^ q is the bolometric luminosity of the sun that according to Cox 
(2000) is 4.74. The equation that regulates the total luminosity of a star 
with it's mass is 

loglO ( 7^ ) = flLM + &LM logio ( YY- ) , (2) 
Lq Mq 

here L is the total luminosity of a star , Lq the sun's luminosity, A4 the star's 

mass , Mq the sun's mass , Olm and 6lm two coefficients that are reported 
in Table 1 for MAIN V , GIANTS III and SUPERGIANTS I ; more details 
can be found in Zaninetti (2005). We remember that the tables of calibration 



Table 1: Table of coefficients derived from the calibrated data (see Table 15.7 
in Cox (2000)) through the least square method 





MAIN, V 


GIANTS, III, 


SUPERGIANTS I 




-0.641 ± 0.01 


-0.792 ± 0.06 


-0.749 ± 0.01 


Tbv[K] 


7360 ± 66 


8527± 257 


8261 ± 67 


Kbc 


42.74 ± 0.01 


44.11 ± 0.06 


42.87 ± 0.01 


Tbc[K] 


31556 ± 66 


36856 ± 257 


31573 ± 67 


olm 


0.062 ± 0.04 


0.32 ± 0.14 


1.29 ± 0.32 




3.43 ± 0.06 


2.79 ± 0.23 


2.43 ± 0.26 



of MK spectral types unify SUPERGIANTS la and SUPERGIANTS lb into 
SUPERGIANTS I, see Table 15.7 in Cox (2000) and Table 3.1 in Bowers 
& Deeming (1984). Prom the theoretical side Padmanabhan (2001) quotes 
3 < ^LM < 5 ; the fit on the calibrated values gives 2.43 < 6lm < 3.43 , see 
Table 1. 
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From a visual inspection of formula (2) is possible to conclude that a 
logarithmic expression for the mass as function of the temperature will allows 
us to continue with formulas easy to deal with. The following form of the 
mass-temperature relationship is therefore chosen 

M T 
^ogwi^-j—) = auT + buT^ogioi—) , (3) 

where T is the star's temperature, Tq the sun's temperature , qmt and 6mt 
two coefficients that are reported in Table 2 when the masses as function of 
the temperature ( e.g. Table 3.1 in Bowers & Deeming (1984)) are processed. 
According to Cox (2000) Tq = 5777 K. Due to the fact that the masses of 



Table 2: Tabic of omt and 6mt with masses given by Table 3.1 in Bowers & 
Deeming (1984). 





main,V 


giants, III 


supergiants,! 










(B-V)> 0.76 


(B-V) < 0.76 




-7.6569 


5.8958 


-3.0497 


4.1993 




2.0102 


-1.4563 


-0.8491 


1.0599 


x' 


28.67 


3.41 


20739 


18.45 



the SUPERGIANTS present a minimum at {B-V)^ 0.7 or T^ 5700 K we 
have divided the analysis in two. Another useful formula is the bolometric 
correction BC 

SC = Mboi-Mv = -^-10 logio T + Xbc , (4) 

where Mboi is the absolute bolometric magnitude. My is the absolute visual 
magnitude, Tbc and Kbc are two parameters that are derived in Table 1 . 
The bolometric correction is always negative but in Allen (1973) the analyt- 
ical formula was erroneously reported as always positive. 

The fifth equation connects the physical variable T with the observed 
color {B — V) , see for example Allen (1973), 

(S-T/)=Xbv + ^ , (5) 

where Kbv and Tbv are two parameters that are derived in Table 1 from 
the least square fitting procedure. Conversely in Section 5 we will explore 
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a series development for {B — V) as given by a Taylor series in the variable 
1/T. Inserting formulas (4) and (3) in (2) we obtain 



My — —2.5 OlM — 2.5 &LM O'MT ~ 
2.5 &LM ^MT ^ioT 

-irBC + lO%ioT+^ + Mb„i,0 . 



(6) 



Inserting equation (5) in (6) the following relationship that regulates My and 
[B — V) in the H-R diagram is obtained 



Up to now the parameters qmt and ^mt are deduced from Table 3.1 
in Bowers & Deeming (1984) and Table 2 reports the merit function 
computed as 



where My represents the calibration value for the three MK classes as given 
by Table 15.7 in Cox (2000). From a visual inspection of the reported 
in Table 2 we deduced that different coefficients of the mass-temperature 
relationship (3) may give better results. We therefore found by a numerical 
analysis the values Omt and buT that minimizes equation (8) when {B — V) 
and My are given by the calibrated values of Table 15.7 in Cox (2000). 
This method to evaluate auT and 6mt is new and allows to compute them 
in absence of other ways to deduce the mass of a star. The absolute visual 
magnitude with the data of Table 3 is 




My = -2.5 olm - 2.5 6lm omt - 
2.5 &LM buT logio{jr= — ) 



(7) 



n 



(8) 



Mv = 31.34- 
3.365 In (7361.0 {{B - V) + 0.6412)"^) + 

4.287(5 - V) 



(9) 
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Table 3: Table of omt and byiT when My'' is given by calibrated data. 





main,V 


giants, III, 


supergiants, I 










(B-V)>0.76 


(B-V)<0.76 




-7.76 


3.41 


3.73 


0.20 




2.06 


-2.68 


-0.64 


0.24 


x' 


11.86 


0.152 


0.068 


0.567 



MAIN SEQUENCE , V when 
- 0.33 <{B-V) < 1.64, 



Mv = -109.6 + 
12.51 In (8528.0 ((5 - V) + 0.7920)"^) 

+4.322(5- F) (10) 
GIANTS , III 0.86 <(B-V)< 1.33 , 



My = -39.74 + 
3.565 In (8261.0 {(B - V) + 0.7491)"^) + 

3.822 (5 -F) (11) 
SUPERGIANTS, I when 
-0.27< (B-V)< 0.76, 



My = -61.26 + 
6.050 In (8261.0 ((5 - V) + 0.7491)"^) + 

3.822 (B - V) (12) 
SUPERGIANTS , I when 
0.76 <(B-V)< 1.80 . 

Is now possible to build the calibrated and theoretical H-R diagram , see 
Figure 1. 
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Figure 1: My against {B — V) for calibrated MK stars (triangles) and theo- 
retical relationship with cmt and 6mt as given by Table 3. 
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2.1 The mass (B-V) relationship 

Is now possible to deduce the numerical relationship that connects the mass 
of the star with the variable {B — V) and the constants omt and 6mt 

>Mrln( (^ J-_^^J (ln(lQ)r' . (13) 

When omt and 6mt are given by Table 3 and the other coefficients are as 
reported in Table 1 the following expression for the mass is obtained: 

log.o(^^) = 7.769 + 

+0.8972 In (7360.9 ({B - V) + 0.6411)"^) (14) 
MAIN SEQUENCE , V when 
-0.33 <(B-V)< 1.64 , 



-1.167 In (8527.5 {{B - V) + 0.792)"^) (15) 
GIANTS , III 0.86 <(B-V)< 1.33 , 

M 

+0.1276 In (8261 {{B - V) + 0.7491)"^) (16) 
SUPERGIANTS , I when 
-0.27 <{B-V)< 0.76 , 



logio(TT^) = 3.73 + 



,M 

-0.2801 In (8261 {{B - V) + 0.7491)"^) (17) 
SUPERGIANTS , I when 
0.76 <{B-V) < 1.80 . 
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against {B - V) for calibrated MK stars : MAIN SE- 



QUENCE , V (triangles) , GIANTS , III (empty stars) and SUPER- 
GIANTS , I (empty circles). The theoretical relationship as given by for- 
mulas (14-17) is reported as a full line. 



Figure 2 reports the logarithm of the mass as function of (B — V) for the three 
classes here considered (points) as well the theoretical relationships given by 
equations (14-17) (full lines). 



2.2 The radius {B - V) relationship 

The radius of a star can be found from the Stefan-Boltzmann law, see for 
example formula (5.123) in Lang (1999) . In our framework the radius is 



log 



10 ^ 



A- 

R&' 



infr©) 

1/2 aiM + 1/2 bLM aMT + 2 -^^^^ + 
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+i/2i..Mi.Mr In [^j^zrvy^) d-d"))"' 

When the coefficients are given by Table 3 and Table 1 the radius is 

logio(^) = -5.793 + 

0.6729 In (7360 {{B - V) + 0.6411)"^) (19) 
MAIN SEQUENCE , V when 
-0.33 <(B-V)< 1.64 , 



logio(;^) = 22.25 

-2.502 In (8527 ((5 - V) + 0.7920)"^) (20) 
GIANTS , /// when 
0.86 <{B-V)< 1.33 , 

logio(^) = 8.417- 

0.7129 In (8261 {{B - V) + 0.7491)"^) (21) 
SUPERGIANTS , I when 
-0.27 <{B-V) < 0.76 , 

log.o(j|) = 12.71 

-1.21 In (8261 {{B - V) + 0.7491)"^) (22) 
SUPERGIANTS , / when 
0.76 <{B-V) < 1.80 . 

Figure 3 reports the radius as function of (B — V) for the three classes (points) 
as well the theoretical relationship given by equations (20-23) (full lines). 
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Figure 3: logio(;^)) against {B - V) for calibrated MK stars : MAIN SE- 
QUENCE , V (triangles) , GIANTS , III (empty stars) and SUPER- 
GIANTS , I (empty circles). The theoretical relationships as given by for- 
mulas (20-23) are reported as a full hne. 
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2.3 The luminosity {B — V) relationship 

The luminosity of a star can be parametrized as 

('"^ '° [(B-V)'-K,y) STm) ■ 
When the coefficients are given by Table 3 and Table 1 the luminosity 

logio(;^) = -26.63 + 

+3.083 In (7360.9 {{B - F) + 0.6411)"^) 
MAIN SEQUENCE , V when 
-0.33 <(B-V)< 1.64 , 



logio(;^) = 29.469 + 

-3.2676 In (8527.59 ((5 - V) + 0.7920)"^) 
GIANTS , III 0.86 <(B-V)< 1.33 , 

logio(-^) = 1.7881 + 

+0.3112 In (8261.19 {{B - V) + 0.7491)'^) 
SUPERGIANTS , I when 
-0.27 <(B-V)< 0.76 , 



logio(-^) = 10.392 + 

-0.682 In (8261.19 ((5 - V) + 0.749)"^) 
SUPERGIANTS , I when 
0.76 <{B-V)< 1.80 , 
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Figure 4: logiol^;^) against (5 - V) for calibrated MK stars : MAIN SE- 
QUENCE , V (triangles) , GIANTS , III (empty stars) and SUPER- 
GIANTS , I (empty circles). The theoretical relationship as given by for- 
mulas (25-27) is reported as a full line. 

Figure 4 reports the luminosity as function of [B — F) for the three classes 
(points) as well as the theoretical relationship given by equations (25-27) 
(full hues). 



3 Absence of calibrated physical parameters 

The already derived framework can be applied to a class of stars in which 
the calibration data of BC , (B — y) , 7W , L versus the temperature are 
absent, for example the white dwarfs. The presence of the fourth edition 
of the Villanova Catalog of Spectroscopically Identified White Dwarfs, see 
McCook & Sion (1999), makes possible to build the H-R diagram of 568 
white dwarfs that have trigonometric parallax. Once the observed absolute 
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magnitude is derived the merit function -y^ is computed as 

n 

X^-Y.^My-Mtf , (28) 

j=i 

where My^ represents the observed value of the absolute magnitude and the 
theoretical absolute magnitude, My , is given by equation (7). The four 
parameters Kby , Tbv , Tbc and Kbc are supplied by theoretical arguments, 
i.e. numerical integration of the fluxes as given by the Planck distribution 
,see Planck (1901). The remaining four unknown parameters omt , &mt , 
Olm and 6lm are supplied by the minimization of equation (28). Table 4 
reports the eight parameters that allow to build the H-R diagram. 

Table 4: Table of the adopted coefficients for white dwarfs. 



Coefficient 


value 


method 


Kby 


-0.4693058 


from the Planck law 


Tbv[K] 


6466.229 


from the Planck law 


Kbc 


42.61225 


from the Planck law 


Tbc[K] 


29154.75 


from the Planck law 


olm 


0.28 


minimum on real data 


^'lm 


2.29 


minimum on real data 


omt 


-7.80 


minimum on real data 


6mt 


1.61 


minimum on real data 



The numerical expressions for the absolute magnitude (equation (7)) , 
radius (equation (18)) , mass (equation (13)) and luminosity (equation (23)) 
are: 

My = 8.199 

+0.3399 In (6466.0 ((5 - V) + 0.4693)"^) + 

4.509(5- F) (29) 
white dwarf, -0.25 <{B-V) < 1.88 , 

logio(7^) = -1-267 

-0.0679 In (6466 ((5 - V) + 0.4693)"^) (30) 
white dwarf, -0.25 < {B - V) < 1.88 , 
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Figure 5: My against (B — V) (H-R diagram) of the fourth edition of the Vil- 
lanova Catalog of Spectroscopically Identified White Dwarfs. The observed 
stars are represented through small points, the theoretical relationship of 
white dwarfs through big full points and the reference relationship given by 
formula (33) is represented little square points. 
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+0.6992 In (6466 {{B - V) + 0.4693)"^) (31) 
white dwarf, -0.25 < {B - V) < 1.88 , 



logio(i—) = -17.58 + 

Lq 

1.601 In (6466 {{B - V) + 0.4693)"^) (32) 
white dwarf, -0.25 <{B -V) < 1.88 . 

Figure 5 reports the observed absolute visual magnitude of the white 
dwarfs as well as the fitting curve; Table 5 reports the minimum , the average 
and the maximum of the three derived physical quantities. 

Table 5: Table of derived physical parameters of the Villanova Catalog of 
Spectroscopically Identified White Dwarfs 



parameter 


min 


average 


maximum 


M/Mq 

R/Rq 
L/Lq 


5.28 10-3 
1.07 10-2 
1.16 10-5 


4.23 10-2 
1.31 10-2 
2.6 10-3 


0.24 

1.56 10-2 
7.88 10-2 



Our results can be compared with the color-magnitude relation as sug- 
gested by McCook & Sion (1999), where the color- magnitude calibration due 
to Dahn et al. (1982) is adopted , 

My = 11.916 {{B -V) + 1)°-^^ - 0.011 
when{B -V) < 0.4 

(33) 

Mv = 11.43 + 7.25 {B-V)- 3.42 {B - Vf 

when 0.4 <{B-V) 

The previous formula (33) is reported in Figure 5 as a line made by little 
squares and Table 6 reports the computed as in formula (28) for our 
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formula (30) and for the reference formula (33). From a visual inspection of 
Table 6 is possible to conclude that our relationship represents a better fit 
of the data in respect to the reference formula. 



Table 6: Table of when the observed data are those of the fourth edition 
of the Villanova Catalog of Spectroscopically Identified White Dwarfs. 

equation 

our formula (30) 710 
reference formula (33) 745 



The three classic white dwarfs that are Procyon B, Sirius B, and 40 
Eridani-B can also be analyzed when {B~V) is given by Wikipcdia (http: / /www.wikipedia.org/). 
The results are reported in Table 7, 8 and 9 where is also possible to visualize 
the data suggested by Wikipedia. 



Table 7: Table of derived physical parameters of 40 Eridani-B where {B — 
y)=0.04 



parameter 


here 


suggested in Wikipedia 


My 


11.59 


11.01 


M/Mq 


6.41 10^2 


0.5 


R/Rq 


1.23 10-2 


2 10-2 


L/Lq 


3.53 10-3 


3.3 lO-''^ 



Table 8: Table of derived physical parameters of Procyon B where {B — 



parameter 


here 


suggested in Wikipedia 


My 


11.43 


13.04 


M/Mq 


7.31 10-2 


0.6 


R/Rq 


1.21 10-2 


2 10-2 


L/Lq 


4.77 10-=^ 


5.5 10-4 
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Table 9: Table of derived physical parameters of Sirius B where {B — V)—- 



0.03 



parameter 


here 


suggested in Wikipedia 


My 


11.32 


11.35 


M/Mq 


8.13 10-2 


0.98 


R/Rq 


1.21 10-2 


0.8 10-2 


L/Lq 


6.08 10-3 


2.4 10-3 



4 Application to the astronomical environ- 
ment 

The stars in the first 10 pc , as observed by Hipparcos (ESA (1997)) , belong 
to the MAIN V and Figure 6 reports the observed stars, the calibration stars 
and the theoretical relationship given by equation (10) as a continuous line. 

Different is the situation in the first 50 pc where both the MAIN V and 
the GIANTS III are present , see Figure 7 where the theoretical relationship 
for GIANTS III is given by equation (11). 

Another astrophysical environment is that of the Hyades cluster, with 
{B — V) and as given in Stern et al. (1995) and available on the VizieR 
Online Data Catalog. The H-R diagram is build in absolute magnitude 
adopting a distance of 45 pc for the Hyades, see Figure 8 , where also the 
theoretical relationship of MAIN V is reported. 

Another interesting open cluster is that of the Pleiades with the data of 
Micela et al. (1999) and available on the VizieR Online Data Catalog. 

Concerning the distance of the Pleiades wc adopted 135 , according 
to Bouy et al. (2006) ; other authors suggest 116 pc as deduced from the 
Hipparcos data , see Mermilliod et al. (1997). The H-R diagram of the 
Pleiades is reported in Figure 9. 

4.1 Distance determination 

The distance of an open cluster can be found through the following algorithm 

1. The absolute magnitude is computed introducing a guess value of the 
distance. 
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Figure 6: My against {B — V) (H-R diagram) in the first 10 pc. The observed 
stars are represented through points, the cahbrated data of MAIN V with 
great triangles and the theoretical relationship of MAIN V with a full line. 
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Figure 7: My against (B — V) (H-R diagram) in the first 50 pc. The observed 
stars are represented through points, the cahbrated data of MAIN V and 
GIANTS III with great triangles , the theoretical relationship of MAIN V 
and GIANTS III with a full line. 
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B-V 



Figure 8: My against {B — V) (H-R diagram) for the Hyades . The observed 
stars are represented through points, the theoretical relationship of MAIN V 
with a full line. 
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Figure 9: My against {B — V) (H-R diagram) for the Pleiades. The observed 
stars are represented through points, the theoretical relationship of MAIN V 
with a full line. 
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2. Only the stars belonging to MAIN V are selected. 



3. The between observed an theoretical absolute magnitude ( see for- 
mula (10)) is computed for different distances. 

4. The distance of the open cluster is that connected with the value that 
minimize the x^. 

In the case of the Hyadcs this method gives a distance of 37.6 pc with an 
accuracy of 16% in respect to the guess value or 19% in respect to 46.3 pc of 
Wallerstein (2000). 

5 Theoretical relationships 

In order to confirm or not the physical basis of formulas (4) and (5) we per- 
formed a Taylor-scries expansion to the second order of the exact equations 
as given by the the Planck distribution for the colors , see Section 5.1 , and 
for the bolometric correction , see Section 5.2. A careful analysis on the 
numerical results applied to the Sun is reported in Section 5.3. 

5.1 Colors versus Temperature 

The brightness of the radiation from a blackbody is 



where c is the light velocity , k the Boltzmann constant , T the equivalent 
brightness temperature and A the considered wavelength , see formula (13) 
in Planck (1901), or formula (275) in Planck (1959), or formula (1.52) in 
Rybicki & Lightman (1985) , or formula (3.52) in Kraus (1986). 
The color-difference , (Ci - C2) , can be expressed as 



where S\ is the sensitivity function in the region specified by the index A , 
i^T is a constant and Ix is the energy flux reaching the earth. We now define 
a sensitivity function for a pseudo-monochromatic system 




(34) 



(Ci - C2) =mi-m2^ K - 2.51ogio 



IS2lxdX 



(35) 



Sx = S{X -Xi) i = U, B, V, R, I 



(36) 
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where 6 denotes the Dirac delta function. In this pseudo-monochromatic 
color system the color- difference is 



{Ci - C2) = K - 2.5log 



(exp( 



10 



he ^ 



XI (exp( 



he • 



1) 



-1) 



(37) 



Table 10: Johnson system 



symbol 


wavelength (A) 


U 


3600 


B 


4400 


V 


5500 


R 


7100 


I 


9700 



The previous expression for the color can be expanded through a Taylor 
series about the point T = 00 or making the change of variable x — ^ , about 
the point x = . When the expansion order is 2 we have 



(Ci — C2)app — 

he {Xi — X2) 5 



In 



A2 



Xiy In(lO) 



A X2X1 kin (10) T 48 A/A/A;2 1n(10)T2 



(38) 



where the index app means approximated. We now continue inserting the 
value of the physical constants as given by CODATA Mohr & Taylor (2005) 
and wavelength of the color as given by Table 15.6 in Cox (2000) and visible in 
Table 10. The wavelength of U,B and V are exactly the same of the multicolor 
photometric system defined by Johnson (1966), conversely R(7000 A) and 
1(9000 A) as given by Johnson (1966) are slightly different from the values 
here used. We now continue parameterizing the color as 



(Ci - C2) 



app 



b d 

a+- + — 



(39) 



Another important step is the calibration of the color on the maximum 
temperature Tcai of the reference tables. For example for MAIN SEQUENCE 
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V at Teal = 42000 , see Table 15.7 in Cox (2000) , {B - V) = -0.3 and 
therefore a constant should be added to formula (38) in order to obtain such 
a value. With these recipes we obtain , for example 

..X . . 3543 17480000 

{B-V)^ -0.4243 + — + (40) 

MAIN SEQUENCE,V when - 0.33 <(B-V)< 1.64. 

The basic parameters b and d for the four colors here considered are re- 
ported in Tabic 10. The parameter a when WHITE DWARF , MAIN SE- 
QUENCE, V, GIANTS, III and SUPERGIANTS, I are considered is reported 
in Table 11 , conversely the Table 12 reports the coefficient b and d that are 
in common to the classes of stars here considered. The WHITE DWARF 



Table 11: Coefficient a 





(B-V) 


(U-B) 


(V-R) 


(R-I) 


WHITE DWARF , T^^^K] =25200 


- 0.1981 


- 1.234 






MAIN SEQUENCE, V , Teai[K] =42000 


- 0.4243 


- 1.297 


- 0.233 


- 0.395 


GIANTS III , Tc„i[K]=5050 


- 0.5271 


- 1.156 


- 0.4294 


- 0.4417 


SUPERGIANTS I , TcaiM =32000 


- 0.3978 


- 1.276 


- 0.2621 


- 0.420 



calibration is made on the values of {U — B) and {B — V) for Sirius B, see 
Wikipedia (http:/ /www.wikipedia.org/). 



Table 12: Coefficients b and d 





(B-V) 


(U-B) 


(V-R) 


(R-I) 


b 


3543 


3936 


3201 


2944 


d 


17480000 


23880000 


12380000 


8636000 



The Taylor expansion agrees very well with the original function and 
Figure 10 reports the difference between the exact function as given by the 
ratio of two exponential and the Taylor expansion in the {B — V) case. 

In order to establish a range of reliability of the polynomial expansion we 
solve the nonlinear equation 

(Ci - C2) - (Ci - C2)„pp = /(T) = -0.4 , (41) 
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Figure 10: Difference between {B — V) , the exact value from the Planck 
distribution , and {B — V)app , approximate value as deduced from the Taylor 
expansion for MAIN SEQUENCE, V 
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Table 13: Range of existence of the Taylor expansion for MAIN SE- 
QUENCE, V 







(B-V) 


(U-B) 


(V-R) 


(R-I) 


T • 


[K] 


4137 


4927 


3413 


2755 


T 

max 


[K] 


42000 


42000 


42000 


42000 



for T. The solutions of the nonlinear equation are reported in Table 13 for 
the four colors here considered. For the critical difference we have chosen the 
value —0.4 that approximately corresponds to 1/10 of range of existence in 
{B — V) . Figure 11 reports the exact and the approximate value of {B — V) 
as well as the calibrated data. 



5.2 Bolometric Correction versus Temperature 

The bolometric correction BC , defined as always negative, is 

BC = M^^,-My , (42) 

where Mboi is the absolute bolometric magnitude and My is the absolute 
visual magnitude. It can be expressed as 



5 "'M "^VfcTTrV Wv' exp(p^)-l 

where Ay is the visual wavelength and Kbc a constant. We now expand with 
a Taylor series about the point T — oo 



15 In (T) _ 5 he _ 5 h^c^ 
Y In (10) ~ 4 kXv ln(10)T ~ 48 k'^Xv^ln{10)T^ 



BCapp 



2,._...^^, ■ (44) 



The constant Kbc can be found with the following procedure. The maximum 
of BCapp is at Tjnax j whcrc the index max stands for maximum 

T = - ^ C451 

max — ^ , X ■ V '-^/ 

D KAv 
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Figure 11: Exact {B — V) as deduced from the Planck distribution , or 

equation (37), traced with a full line. Approximate {B — V) as deduced 
from the Taylor expansion, or equation (38), traced with a dashed line. The 
cahbrated data for MAIN SEQUENCE V are extracted from Table 15.7 
in Cox (2000) and are represented through empty stars. 
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Figure 12: Exact BC as deduced from the Planck distribution , or equa- 
tion (43), traced with a fuU hne. Approximate BC as deduced from the Tay- 
lor expansion, or equation (46), traced with a dashed line. The calibrated 
data for MAIN SEQUENCE V are extracted from Table 15.7 in Cox (2000) 
and are represented through empty stars. 

Given the fact that the observed maximum in the BC is -0.09 at 7300 K in 
the case of MAIN SEQUENCE V we easily compute Kbc and the following 
approximate result is obtained 

= 31.41 -3.257 in (D- if" . (46) 

Figure 12 reports the exact and the approximate value of BC as well as the 
cahbrated data. 

The Taylor expansion agrees very well with the original function and Fig- 
ure 13 reports the difference between exact function as given by equation (43) 
and the Taylor expansion as given by equation (46). 

In order to establish a range of reliability of the polynomial expansion we 
solve the nonlinear equation for T 

BC-BCapp = -OA . (47) 



29 



o 




Figure 13: Difference between (B — V) , the exact value from tlie Planck 
distribution , and {B — V)app , approximate value as deduced from the Taylor 
expansion for MAIN SEQUENCE, V 
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Table 14: (B-V) of the Sun , T=5777 K 



meaning 


(B-V) 


calibration, Cox (2000) 


0.65 


here, Taylor expansion 


0.711 


here , Planck formula 


0.57 


least square method, Zaninetti (2005) 


0.633 


Allen (1973) 


0.66 


Sekiguchi & Fukugita (2000) 


0.627 


Johnson (1966) 


0.63 



Table 15: (R-I) of the Sun , T=5777 K 



meaning (R-I) 



cahbration. Cox (2000) 


0.34 


here, Taylor expansion 
here , Planck formula 


0.37 
0.34 



The solution of the previous nonlinear equation allows to state that the 
bolometric correction as derived from a Taylor expansion for MAIN SE- 
QUENCE, V is rehable in the range 4074 K <T < 42000 K. 

5.3 The Sun as a blackbody radiator 

The framework previously derived allows to compare our formulas with one 
specific star of spectral type G2V with T=5777 K named Sun. In order 
to make such a comparison wc reported in Table 14 the various value of 
[B — V) as reported by different methods as well as in Table 15 the value 
of the infrared color (R-I). From a careful examination of the two tables we 
conclude that our model works more properly in the far-infrared window in 
respect to the optical one. 
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6 Conclusions 



New formulas A new analytical approach based on five basic equations 

allows to connect the color {B — V) of the stars with the absolute visual 
magnitude, the mass, the radius and the himinosity.The suggested method is 
based on eight parameters that can be precisely derived from the calibration 
tables; this is the case of MAIN V , GIANTS III and SUPERGIANTS I. In 
absence of cahbration tables the eight parameters can be derived mixing four 
theoretical parameters extracted from Planck distribution with four parame- 
ters that can be found minimizing the connected with the observed visual 
magnitude ; this is the case of white dwarfs. In the case of white dwarfs the 
mass-luminosity relationship , see Table 4, is 

logio(^) = 0.28 + 2.29 logio(^) , (48) 
white dwarf, Q.QQ5Mq < M < /.gAMq . 

Applications The applications of the new formulas to the open clusters 
such as Hyades and Pleiades allows to speak of universal laws for the star's 
main parameters. In absence of accurate methods to deduce the distance of 
an open cluster an approximate evaluation can be done. 

Theoretical beises The reliabihty of an expansion at the second order 
of the colors and bolometric correction for stars as derived from the Planck 
distribution is carefully explored and the range of existence in temperature 
of the expansion is determined. 

Inverse function In this paper we have chosen a simple hyperbohc 
behavior for {B — V) as function of the temperature as given by formula (5). 
This function can be easily inverted in order to obtain T as function of (B—V) 
(MAIN SEQUENCE, V) 

T=-, K (49) 

{B-V) + 0.641 ^ ' 

MAIN SEQUENCE, V when 4137 < T[K] < 42000 
or when - 0.33 < {B - V) < 1.45 . 

When conversely a more complex behavior is chosen , for example a two 
degree polynomial expansion in ^ as given by formula (38) , the inverse 
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formula that gives T as function of {B - V) (MAIN SEQUENCE, V) is more 
complicated than formula (49), 



The mathematical treatment that allows to deduce the coefficients of the 
scries reversion can be found in Morse & Feshbach (1953); Dwight (1961); 
Abramowitz & Stegun (1965). 
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